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ABSTRACT 

We find a linear instability of non-axisymmetric Rossby waves in a thin non- 
magnetized Keplerian disk when there is a local maximum in the radial profile of a 
key function £,(r) = J-(r)S 2 / r (r), where J~~ x = z • (V x v)/S is the potential vorticity, 
S = P/YF is the entropy, S is the surface mass density, P is the vertically integrated 
pressure, and T is the adiabatic index. We consider in detail the special case where 
there is a local maximum in the disk entropy profile S(r). This maximum acts to 
trap the waves in its vicinity if its height to width ratio max(5)/Ar is larger than a 
threshold value. The pressure gradient derived from this entropy variation provides the 
restoring force for the wave growth. We show that the trapped waves act to transport 
angular momentum outward. A plausible way to produce an entropy variation is 
when an accretion disk is starting from negligible mass and temperature, therefore 
negligible entropy. As mass accumulates by either tidal torquing, magnetic torquing, 
or Roche-lobe overflow, confinement of heat will lead to an entropy maximum at the 
outer boundary of the disk. Possible nonlinear developments from this instability 
include the formation of Rossby vortices and the formation of spiral shocks. What 
remains to be determined from hydrodynamic simulations is whether or not Rossby 
wave packets (or vortices) "hold together" as they propagate radially inward. 



Subject headings: Accretion Disks — Hydrodynamics — Instabilities — Waves 



1. Introduction 

The central problem of accretion disk theory is how angular momentum is transported 
outward so that matter can accrete onto the central gravitating object. Ever since Shakura 
and Sunyaev (1973) and Pringle (1981) outlined the basic structure of accretion disk (see also 
Frank, King, fe Raine 1985| ), the origin of enhanced (or a) viscosity has been the subject of 
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hundreds of research papers. Classically, high Reynolds number turbulence has been invoked to 
explain the enhanced viscosity. The driving mechanism(s) of the turbulence in high-temperature 
(for example, Active Galactic Nuclei) disk systems has been a long standing problem because 
of the apparent absence of linear instabilities in non-magnetized Keplerian disks (see review by 
Papaloizou fc Lin 1995| ). For this reason, magnetic fields and associated instabilities have been 



suggested as the origin of the enhanced viscosity in disks ( Balbus fc Hawley 1998 ). The widely 



cited magneto-rotational instability ( |Velikhov 1959fc phandrasekhar 1960 ) occurs only when the B 



field is weak. Also, it has been studied only for small patches of a disk using unphysical boundary 
conditions on the magnetic field. The nature of the global angular momentum transport remains 
an open problem. 

This paper reconsiders the question of the stability of non-magnetized Keplerian disks. In 
particular, we find a linear instability of non-axisymmetric perturbations for conditions where the 
disk quantities, such as surface density and entropy have steep radial gradients. The conditions 
we consider are in general nonbarotropic which distinguish our work from that of Papaloizou and 
Pringle (1984, 1985; Goldreich, Goodman, & Narayan 1986; Narayan, Goldreich, & Goodman 
1987). Also, in contrast with the work of Papaloizou and Pringle, the modes we consider are 
trapped at least initially in a narrow range of radii and therefore do not depend on reflections 
from inner and outer radii of the disk (or tori). 

The instability we find will give rise to Rossby vortices in the nonlinear limit. These vortices 
are well-known in thin, two dimensional planetary atmospheres, for example, Jupiter's "Red 
Spot." The persistence of such vortices would be crucial for the hydrodynamic transport of angular 
momentum in accretion disks. Here, we present our findings of initial conditions which cause 
Rossby waves to grow, and defer the nonlinear behavior of the vortices to future publications. 

We emphasize that the astrophysical relevance of this instability is the assumption that in 
the circumstances where accretion disks are important, mass at low temperature and entropy 
accumulates at some external radius to a thickness such that a condition of a local maximum 
in entropy is created because of confinement of heat. Then a wave of non- linear Rossby vortex 
excitation carries the mass, and entropy maximum inwards, exciting further vortices that transport 
the angular momentum outwards. 

In §2, we consider the general stability of a 2D nonbarotropic disk and point out the role of 
a key function C{r) in determining the stability. We then study a case with a significant 'bump' 
in the radial distribution of entropy (or temperature) and calculate its growth rate and dispersion 
relation. In §3, we discuss the possible nonlinear evolution of Rossby waves, especially their role 
in angular momentum transport. 
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2. Theory 

We consider the stability of a thin (2D), non-magnetized Keplerian accretion disk. In treating 
the disk as 2D we assume that the vertical thickness 2h is <C r, perturbation wavelengths are 
much longer than 2h. Using an inertial cylindrical (r, 0, z) coordinate system, we have a surface 
mass density S(r) = J h h dzp(r, z), and a vertically integrated pressure P(r) = J^ h dzp(r, z), We 
consider the equilibrium disk to be steady (d/dt = 0) and axisymmetric (d/d(f) = 0), with the flow 
velocity v = v^(r)^> (that is, the accretion velocity v r is neglected), the gravitational potential 
<J>(r) ps —GM/r (the self-gravity is neglected in this study), and the radial force equilibrium, 
Sv^/r = dP/dr + Y>d$/dr, where M is the mass of the central object. The vertical hydrostatic 
equilibrium gives h ps {c s /v^)r, where c s = (rP/E) 1 / 2 is the sound speed and T is the adiabatic 
index. 



2.1. Basic Equations for a Nonbarotropic Disk 

The perturbations considered are in the plane of the disk so that the perturbed surface mass 
density is E = E + <5S(r, (f>, t); the perturbed vertically integrated pressure is P = P + 5P(r, 4>, t); 
and the perturbed flow velocity is v = v + 5v(r,(f),t) with 5v = (<5tv, <k^, 0). The equations for 
the perturbed disk are 

^ + SV-v = 0, (la) 
Dv 1 

_ = _VF-V 4 », (16) 



D 
Dt 



P 



sr , (lc) 

where D/Dt = d/dt + v • V and we refer to S = P/TF as the entropy of the disk matter. Equation 
(lc) corresponds to the isentropic behavior of the disk matter. We use equation (lc) for simplicity 
rather than the rigorous relation D(p/ p y )/Dt = (see Papaloizou and Lin 1995). 

We consider perturbations oc /(r)exp(im0 — iujt), where m is the azimuthal mode number 
and uj the angular frequency. From equation (la), we have 

iAco <$E = V • (E <5v) , (2) 

where Aw(r) = lo — mQ(r) and 17 = v^/r ps (GM/r 3 )i From equation (lb) we have 

iAu8v r + 2VL5 V<t> = -— - ^ -fo , (3°) 
iAuSvtj, - -j^5v r = ik^— . (3b) 
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Here, k = [r ^ d(r A £l 2 ) / dr] 2 « f2 is the radial epicyclic frequency and fc^ = m/r is the azimuthal 
wavenumber. From equation (lc), we have 



5P = c 2 <5£ - ^-5v r , 



which is the perturbed equation of state and 



L a = r 



dr 



In 



P 

ST 



(4) 



(5) 



is the length scale of the entropy variation in the disk. The nonbarotropic nature of the disk 
renders L s to be finite whereas for a barotropic disk (i.e., P is a function of S only), L s — > 00. 

To further simplify the momentum equation, we can write the right-hand-side of equation 
(3a) as (6P/E)' + (oP/£ 2 )X/ - (o£/£ 2 )P', where the prime denotes d/dr. Next, equation (4) can 
be used to express this as (<5P/£)' — 5P/(EL S ) — i5v r c 2 / '(A00L S L P ) , where 



I ln(P) 



(6) 



is the length-scale of the pressure variation. For definiteness we assume h < (\L S \,\L P \ 
Using * = <5P/£ as our key variable, we can now rewrite equations (3) as 



i\ Alo + 



AujL s L p 



6~iv + 2VL5vs = — j— 



dr 



Solving equations (7) for {5v r ,8v^) gives 



Alo 



— ( *' - — ] - 2^* 



(7a) 
(76) 

(8a) 



+ 



AujUL s L p 



^ + 



+ 



where 



£ 



K 2 _ A(J 2 _ c 2/( LsLp ) 



(86) 
(9) 



For 'corotation modes' where |Acj| 2 is small compared with k 2 and c 2 /\L s L p \, we have 
T Sf]//c 2 = T,/2lo z . Here, w 2 = z ■ (V x v) is the vorticity. The ratio w z /£ « 1/(2JF) is 
commonly referred to as the potential vorticity (see for example Drazin 1978). 
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We next use the continuity equation (2) and equations (4) and (8) to obtain 



IT 



n 



+ 



T 1 



, £* 

= — H r — * 

Aw 



H h 

AwL, 



(10) 



Equation (10) allows the determination of lu and identification of mode structure for general 
disk flows, both stable or unstable. If the flow is barotropic (as has been studied extensively in 
the literature), all the terms on the lower line of this equation vanish, \L S \ — > oo. In this paper, we 
explore additional effects which arise when \L S \ is finite. 



2.2. Instability Condition 



A useful quadratic form can be obtained by multiplying equation (10) by ** and integrating 
over the disk. Assuming ^^'rF/Q. — > for r — > 0, oo, we obtain 

~Jd 2 r^ (\V\ 2 + klW)=Jd 2 r ||*| 2 + 



/ 



d 2 r 



T 1 / rT x ' 

+ 



ML 2 r \ nL 



■J 



2 d 2 r 



Auj AluLc 



I*! 2 + 



\w + 



d 2 r 



Aoj 2 £ILcL 



1*1 



(11) 



Note that only the last two integral terms on the right-hand side of this equation involve the 
complex variable Aw(r). 

Equation (10) can be solved by several methods (e.g., by a relaxation method or by shooting), 
but we defer the full discussion of its solutions to a future publication. In order to bring out 
the essential aspects of equation (10), we make simplifying assumptions which allow an analytic 
solution. First, when \Au>\ is a finite fraction of Q and c s /v^ <C 1, we find that the integral 
involving Aoo 2 is small compared to the integral with Alo, so that the last term in equation (11) is 
negligible. Furthermore, for the 'corotation modes' where |Aa;| 2 /fi 2 <C 1, and for conditions where 
the pressure variation with r is such that k 2 /Q 2 is positive, the function T is positive and real. 
It is physically possible to have k 2 positive and yet have a finite \L S \ because the pressure and 
entropy variations need not to be the same. In fact, this is precisely a feature of nonbarotropic 
flow (see next section). Thus, the terms containing Auj in (11) can be combined to give 

ksF 



(12) 
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Putting equation (12) back into (11) and taking the imaginary part of (11), we find that a 
solution is possible only for conditions where the imaginary part of (12) vanishes (see also Lovelace 
& Hohlfeld 1978; hereafter LH). Thus, we find that an instability is possible only if In 

vanishes at some r, in other words, the key function C(r) = J-(r)S 2 / r (r) has a maximum or 
minimum. 

This newly defined function C(r) warrants further discussion. We believe that this expression 
is new and more general than those given in previous studies. With the above assumptions, we 
can write C(r) as 

C(r) = TS 2 / T « S - f /F oc % £2/r-i T 2/r (13) 

where f2 2 (r) and K 2 (r) are of course constrained by the radial force balance in the equilibrium. 
An extremum of C(r) necessary for instability could come from several sources: (a) In the 



limit of incompressible and isothermal flow (cf. Staley fc Gall 197E for an application to flow in a 



tornado); an extremum in the radial profile of £l(r) is required to have unstable vortices, (b) In 
the studies of LH and Sellwood & Kahn (1991), a cold disk was assumed and S(r) was effectively 
taken as a constant. In this case C(r) is simply T = S/O, which must have a maximum or a 
minimum for instability. Similarly, in the cases considered by Papaloizou & Pringle (1984, 1985), 
where a polytropic (barotropic) P = -P(S) equation of state is assumed, S(r) is again a constant 
throughout the disk so that a maximum or minimum of C(r) oc J 7 is needed for instability, (c) In 
the present study, we treat the more general case where S(r) varies across the disk. An extremum 
in C{r) can result from having an extremum in S{r) directly (that is, a local "bump" in the 
temperature profile of the disk), or by having a large enough gradient in S(r) which causes an 
extremum in the epicyclic frequency k 2 owing to the radial force balance. Furthermore, we find 
below that instability occurs when a threshold in the variation of C(r) is exceeded. In view of 
equation (13), a variation of T(r) is more important in giving a change in C(r) than a comparable 
fractional variation of £(r). 



2.3. Thermodynamic Driving of Vorticity 



It is straightforward to obtain the vorticity equation (cf. Pedlosky 1987 ; Spiegel 1993|) by 
taking the curl of equation (lb) and combining with (la), 



D 
Dt 



VS x VP 

S3 



(14) 



where again, uj z = z ■ V x v is the vorticity. The right-hand side of this equation represents the 
thermodynamic driving of the vorticity. For a barotropic flow, P = P(S), 



LO- 



(15) 
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This is the case for most previous studies (cf. Papaloizou & Pringle 1984; Papaloizou &j 
Pringle 1985| ; |Goldreich, Goodman, fc Narayan 1986| ; |Narayan, Goldreich, fc Goodman 1987] ). 



Equation (15) says that each fluid element conserves its specific vorticity. In contrast, the term 
VS x VP oc VT x VS destroys this conservation and allows pressure forces to "produce" vorticity 
in the flow. In the non-barotropic case surfaces (or lines) of constant surface density and surfaces 
of constant specific entropy do not coincide. 



2.4. Dispersion Relation 

For the 'corotation modes' (| Aoj| 2 <C £l 2 ) and J-(r) ~ S/O, an approximate dispersion 
relation follows from equation (10) by taking the r— dependence ^ oc exp(i/c r r) with (k r r) 2 ^> 1, 
(k r L s ) 2 S> 1, and (k r L p ) 2 ^ 1. The dominant terms in (10) are: 

n 2 c 2 { ' + Aco 2 n 2 L s L p ' 1 j 

or 



Aw - k ^ /Q 



1 + k 2 /i 2 



(ln£)'±J[(m£)'P l+k " ; '" 



L s L p 



(17) 



where k 2 = & 2 + /c? and /i = c s /Q. The dispersion relation (17) is the disk analogue of the Rossby 
wave dispersion relation (see for example Brekhovskikh & Goncharov 1993). Note in particular 
that the group velocity of the waves in the corotating reference frame is in the direction opposite 
to that of the disk rotation. Instability is possible if [(ln£)'] 2 < (1 + k 2 /i 2 )/(L s L p ), which requires 



L s L p > 0. 



From equation (17), the maximum growth rate occurs for (ln£)' = 0, and it is 

max(wj) l/c^lc 2 /^ 2 



n (l + k 2 h 2 )^(L s L 




where vjc is the Keplerian velocity and coi = Im(w) is the growth rate. The maximum growth 
rate is < f2, for a thin disk where h/r ~ c s jv$ <C 1 and for L s L p > h 2 . Note in particular the 
geometrical optics relation for a wavepacket, 

— - -— ~ -k O ^ --— -0 fl9l 

dt ~ dr ~ 2 ^ ' dt ~ dfi~ ' { ) 

so that fc r ~ ifc^Ot+const., and fc^ = const. Thus, a wave packet initially with |fc r /i| <C 1 is 
rapidly sheared by the differential rotation. In a time of order tg^^ = (fA^/tfi) -1 the wave 
evolves into having a very short radial wavelength, \k r h\ > 1, trailing spiral wave. This is discussed 



-8- 



in detail by Goldreich and Lynden-Bell (1965). The maximum amplification a wave can have 
is A = exp(/ dt Wj) ~ exp[max.(u>i)t s h ear ] « exp[(2/3)(c s /t>;f )(r 2 /L s L p ) 2]. The amplification is 
A = O(l) and is probably not significant for L s L p > h? . 

2.5. Trapped Modes 

Here, we consider conditions where C ~ (E/£l)S 2 / r has a maximum or minimum as a function 
of r at ro of width Ar <C ro- The case where F(r) = £/f2 has an extremum was treated earlier 
by LH. Therefore, we consider the case where T ~ const and where S(r) has an extremum. The 
extremum of S(r) is thus an extremum of T(r). As mentioned earlier, a variation of T(r) is more 
important in giving a change in £(r) than the variation of X(r). For a maximum in S*(r), we find 
using equations (17) and (19) that the radial group velocity of the waves v gr = dtu/dk r is directed 
towards ro, whereas for a minimum the waves propagate away from ro. Thus, a maximum of S(r) 
can give trapping of the waves in the vicinity of ro. This trapping obviates the need for reflecting 
disk boundaries considered by Goldreich et al. (1986). 

For definiteness we take 

S(r) = S! + (S 2 - Si)exp[-(r - r ) 2 /(2Ar 2 )] , (20) 

where S^S^ro, and Ar are constants. For simplicity we consider Ar 2 <C r 2 . In this limit, the 
real and imaginary parts of equation (10) for ^ = 5P/Y, = ^ r + i^i can be written as 

< = (K\ VP, - {K 2 )i , 

= (K\ + {K 2 )i * r . (21) 

Here, we have omitted the detailed expressions for (K 2 ) r and (K 2 )i, both of which are real 
functions of L s , Cs/v^, Au/Q, and m (see figures). It can be shown that (K 2 ) r and (K 2 )i possess 
the following properties (assuming uj r /[m£l(ro)] = 1): (K 2 )i is an odd function of r — ro and is 
— > when (r — ro) 2 S> Ar 2 ; but (K 2 ) r is an even function of r — ro. (K 2 ) r is < for r — > ro, 
which corresponds to a classically allowed region of motion in a potential well, and is > for 
(r — ro) 2 ~> Ar 2 , which corresponds to a forbidden region. It is this potential well give rise to 
"bound states" which determines the growth rates u>i. 

The existence of wave growth (uji > 0) and its dependence on S2/S1 can be obtained by first 
considering the potential well for u>i — > 0. The first appearance of a bound state in the well is 
given by the Bohr-Somerfeld quantization condition, 

/(S2/S1) = fdr(-K 2 )? = 7T . (22) 

This condition corresponds to a critical value of S2/S1 so that there can be a bound state with 
Wj/fi > (which is the middle curve in Figure 1). This is physically the threshold required for 
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the growth of Rossby waves. Again, we emphasize that the existence of wave growth depends 
critically on the nonbarotropic nature of disk, because the main contribution to this potential 
well comes from Ak&J-ty /(AloL s ), which is always nonpositive since lZe(l/Atu) oc r — tq whereas 
1/L S oc ro — r. This term vanishes when L s — > oo. 

Figure 2 shows sample results for (K 2 ) r ^ and ^> r ^ by solving equations (21) using a shooting 
method. We omit discussing the detailed numerical method, but only point out that one can make 
use of the symmetry of the potential well as well as the boundary conditions at r = ro and r — > oo. 
The key physical variables are S2/S1 and oJi/Q. We find u>i to have the approximate dependence 
tOi/0, « fei[(S , 2/<S'i) — k2\ 2 ^ for S2/S1 > k2, where ki is the mentioned critical value and k\ is a 
constant. 

Figure 3 shows the vortical nature of the velocity field for the same unstable case as Figure 2. 



2.6. Angular Momentum Transport 

To calculate the flux of angular momentum, we first note that, with the complex factors 
exp(im(j) — iut) removed, 5v r = 5v® + i5v^ , 5v$ = 5v® + i5v^ and <5£ = 5T, E + i5YP , where the 
E, O superscripts denote even or odd real functions of r — ro- Such properties are obtained by 
expressing them using equation (1) and ^ = + i^f ■ Thus across a circle of radius ro, we have 

c 2tt 



rZTT 

Fj(r ) = dtprl fte(£):Ke((5tv)fte(ty 
Jo 



nr 2 S 6vF(r ) 5v*(r ) , 

^ + 7^^- *r " TT^-^i ■ (23) 



V 2Vtm Vro V V r„ 



We find in general that m^^(ro) < 0. Further, we find that each of the two factors in parentheses 
in (23) are positive. Therefore, the trapped wave causes the outward flow of angular momentum 
across r = ro- Furthermore, one notes that the positive transport is from a "nonlinear" term 
(SvrSVff,) since the first order term vanishes. 



3. Discussion 

The astrophysical importance of this Rossby wave instability is the implication that in the 
nonlinear limit, multiple Rossby vortices will form after these waves break and coalesce in a 
few revolutions. Hence, they provide a plausible explanation of the hydrodynamic transport of 
angular momentum in thin, Keplerian accretion disks. The critical condition for the inception of 
this hydrodynamic transport, as derived in this linear case, is a minimum contrast (~ x3) in the 
radial distribution of entropy. Thus, it is important that heat (and thus entropy) be preserved 
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against cooling by radiation in at least tens of revolutions. Roughly speaking, recognizing 
the large dynamic range and uncertainty of opacities, this implies a minimum mass 'thickness' 
E ~ 100— 1000 g cm -2 . An important implication from this minimum mass 'thickness' requirement 
is the prediction on the turn-on (active state) and turn-off (quiescent state) conditions for the 
angular momentum transport by the Rossby vortices. Specifically, we could envisage the following 
scenario: as mass accumulates in the disk by accretion, the resulting trapping of heat allows the 
temperature to rise by the energy input of accretion and an episode of Rossby vortex accretion 
ensues. Depending upon the heating, opacity conditions, and mass accretion rates, there will 
be a radius, a thickness, and therefore a mass where the Rossby instability initiates. This mass 
and mass accretion rate determine the event luminosity and episode rate. Finally we note that a 
given episode terminates when its mass is accreted presumably rapidly before additional mass is 
supplied. This understanding of the turn-on and turn-off conditions for Rossby vortex transport 
may prove relevant for several astrophysical systems, such as the planet formation after the solar 
nebula has cooled (Cameron 1962) and outbursts from CVs and stellar mass black holes (Vishniac 

1997) . 

Although a localized high entropy region (S2/S1 ~ 3) may appear to be a large deviation in 
S(r), the resultant pressure gradient force is still very small compared to the gravitational force. 
In the case represented by all the figures, we have used an unperturbed c s jv^ = 0.045, which gives 
Cs/v^ ~ 0.078 at the peak of perturbed region (S2/S1 ~ 3). But this deviation in S(r) is sufficient 
to drive the growth of non-axisymmetric perturbations. A possible mechanism of formation of 
an entropy bump is from the interaction between the matter from Roche lobe overflow and the 
accretion disk. A spiral shock can form in this situation (Sawada et al. 1986). 

Bracco et al. (1998) have studied barotropic disks with initial perturbations in the potential 
vorticity (u^/E) profile. They showed that vortices can have a rather long life time. Another 
possible nonlinear effect is the formation of spiral shocks which heat the disk and thereby maintain 
the bump against radiative cooling. In either case, requiring a finite entropy bump for the 
instability to grow implies a maximum cooling rate above which a bump which develops in the 
disk quickly decays away and such an enhanced viscosity model turns off. What remains to 
be determined is whether the bump "holds together" as it propagates radially inward. If these 
Rossby vortices are indeed playing a major role in angular momentum transport, then they are 
fundamentally different from the usual a viscosity prescription which is based on homogeneous 
turbulence. In contrast, Rossby vortices are large scale, coherent structures. 

Preliminary numerical simulations using a high resolution PPM code (Nelson et al. 1998) 
appear to confirm the existence of a threshold on S2/S1 for wave growth and the formation of 
anti-cyclonic vortices. Numerical studies are now underway to follow the nonlinear evolution of 
the vortices in nonbarotropic disks and to study the associated transport properties (Nelson et al. 

1998) . 

We thank the referee's insistence in clarifying many of issues which greatly improved this 
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Fig. 1. — Sample profiles of (K 2 ) r in units of 1/Vq are shown for different values of S2/S1. For this 
plot m = 5, Ar/ro = 0.05, T = 5/3, unperturbed c s jv^ = 0.045, L s = L p , and cjj/n = 10~ 3 . The 
spike at r/ro = 1 is due to the finite value of cjj. For the intermediate curve, the Bohr-Somerfeld 
condition (10) is satisfied to a good approximation. For other values of m between 1 and 10, the 
critical value for the bound state is S2/S1 ~ 2.7 + 0.83(m/10)2 . 
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Fig. 2. — Sample profiles of (K 2 ) r ^ and * r> i, obtained by a shooting method, are shown for 
LOi/n = 0.1, m = 5, and Ar/r = 0.05. For this case, S 2 /Si « 3.744 and tfj(r ) « -6.543. 
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Fig. 3. — Velocity field of the unstable perturbation shown in Figure 2, obtained from equations 
(8) and (21). 



